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ABSTRACT 

A high level of complex structure, or "granularity" , has been observed in the distribution of infrared- 
obscuring material towards the Galactic Centre (GC), with a characteristic scale of 5" — 15", correspond- 
ing to 0.2 — 0.6 pc at a GC distance of 8.5 kpc. This structure has been observed in ISAAC images which 
have a resolution of ~ 0.6", significantly higher than that of previous studies of the GC. 

We have discovered granularity throughout the GC survey region, which covers an area of 1.6° x 0.8° 
in longitude and latitude respectively (300 pc x 120 pc at 8.5 kpc) centred on Sgr A*. This granularity is 
variable over the whole region, with some areas exhibiting highly structured extinction in one or more 
wavebands and other areas displaying no structure and a uniform stellar distribution in all wavebands. 
The granularity does not appear to correspond to longitude, latitude or radial distance from Sgr A*. We 
find that regions exhibiting high granularity are strongly associated with high stellar reddening. 

Subject headings: dust, extinction — Galaxy: center — infrared: stars — ISM: structure 
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1. INTRODUCTION 

Observations towards the Galactic Centre (GC) are ex- 
tremely difficult because of high levels of extinction from 
intervening material in the Galactic Plane. The gas and 
dust which form the inter-stellar medium (ISM), if dis- 
tributed homogeneously along our line of sight, would 
produce ^ 100 magnitudes in visual extinction (Laun- 
hardt et al 2002), in stark contradiction to observations 
(e.g. Rieke & Lebofsky 1985; Catchpole et al 1990; Blum 
et al 1996; Dutra et al 2003) which indicate the true av- 
erage extinction is Ay < 30. The explanation inferred 
for this, supported by numerous observations at differ- 
ent wavelengths (Catchpole et al 1990; Lis & Carlstrom 
1994; Dutra et al 2003) is that the ISM towards the GC 
is distributed non-uniformly. Hitherto, no extinction map 
towards the GC has had an angular resolution < 1'. 

We describe the results derived from new, high- 
resolution near-IR images of 26 fields within the nuclear 
bulge. These VLT-ISAAC fields (Bandyopadhyay et al. 
2005), are 2.5 x 2.5 arcmin^ in size, with a plate scale of 
0.1484" per pixel. The fields are distributed throughout 
an area 1.6° x 0.8° (/,6) centred on Sgr A*, avoiding ar- 
eas of known high star formation. We obtained a limiting 
magnitude of J = 23 {S/N = 5), H = 21 and = 20 
{S/N = 10), taken on nights with seeing < 0.6". 

The purpose of this Letter is to report the discovery and 
quantify the size scales of structure in the dust distribu- 
tion unresolvable in previous studies. A follow-up paper 
(Gosling et al. in prep.) presents extinction maps and 
a detailed analysis of comparisons with previous, lower- 
resolution maps such as that of 2MASS and DENIS. 

2. GRANULARITY 

Examination of the VLT-ISAAC fields of the GC reveals 
that the degree of reddening seen in colour-colour (C-C) 
and colour-magnitude (C-M) diagrams of the field popu- 
lations varies widely (see Fig. lb and e). This property is 
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not evident in all fields. C-C plots for some of the fields 
(such as Fig. lb) show two loci of stars, one for foreground 
stars, and another with reddening consistent with tradi- 
tional values for the GC. C-C plots of other fields (such as 
Fig. le) show the locus of the non-local stellar population 
to be greatly extended along the reddening vector, and 
containing stars with reddening considerably in excess of 
the traditional GC reddening values. 

Visual inspection shows that fields where the reddening 
appears to vary substantially within the field also display 
structure in the stellar distribution on sub-arcmin scales. 
This structure takes the form of filaments and clumpy re- 
gions within which the stellar density is dramatically re- 
duced compared with the field average. We present evi- 
dence in §3.3 that these structures are related to the in- 
creased and varying reddening in these fields. The fine 
scale of the structures we have discovered were inaccessi- 
ble to previous lower-resolution surveys such as 2MASS. 

To quantify "granularity", we used statistical analyses 
to measure the structure of the observed stellar distribu- 
tion as a function of wavelength. Wall & Jenkins (2003) 
present techniques for statistical analysis of 2-dimensional 
distributions of points. Their method is used to detect 
structure in galaxy distributions over the whole sky. To 
apply this test it is assumed that the underlying stellar 
distribution is random, and that the cause of the struc- 
ture arises from the intervening extinction. The i^Tg-band 
images will be less affected by extinction than the J- and 
i?-bands. If the underlying distribution of stars is ran- 
dom, then the Ks images should show a measure closer to 
random than equivalent images in J and H . 

The test divides each field into a grid of cells, whose 
sizes were systematically varied from 150" to 1" to quan- 
tify the size scale of the structures. The mean star density 
per square arcsec was determined by dividing the number 
of stars in the entire field by the field area. The deviation 
of the number of stars observed in a grid cell was com- 
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Fig. 1. — a) Kb image of a field with no apparent structure in the stellar distribution, b) C-C diagram of the stars in the field shown in (a). 
There are two main loci of stars: the local population to the bottom-left and the GC population in the centre, c) The measure of granularity 
for the field shown in (a). In all three bands it does not deviate from zero, as expected for a random distribution, d) Kb image of a field with 
obvious regions of low stellar density compared to the field average, e) C-C diagram for the field shown in (d). Note, compared to the C-C 
diagram in (b), the locus for the GC stars is extended to high reddening, f) The measure of granularity for the field shown in (d) shows that 
there is measurable structure in all three bands. 



pared to the mean number per cell. The deviations from 
the mean were summed and the variance /i2 of all of the 
cells was calculated. Because the mean number of stars 
per unit area is a constant, and it is the size of a cell that 
is altering, the variance is normalised to the size of the 
grid cell, as expressed in the formula below. 



y 



(1) 



where y is the measure of the deviation from randomness 
for the whole field, i.e. granularity, and n is the expected 



number of stars in one grid cell. 

The number of stars ranges from 1211 to 5592 in J im- 
ages, 1792 to 6515 in H and 3005 to 9077 in K^. To act as 
comparison data-sets, we produced 400 test fields with an 
equivalent spatial distribution to the observed fields using 
a computer random mmrber generator to position points in 
a grid, 100 each with 1000, 3000, 5000 and 7000 points. We 
applied the same statistical test to these simulated fields 
as to the GC fields. Comparison of the simulated and ac- 
tual data allowed us to determine whether the degree of 
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granularity measured in a GC field, y, was simply a result 
of low stellar density and as such a statistical effect, or the 
result of a non-random apparent stellar distribution. The 
maximum values of y for the random fields were 1.8, 0.08, 
0.05 and 0.04 for the fields with 1000, 3000, 5000 and 7000 
points respectively, at a scale of 2" (Fig. 2). 

3. RESULTS 

3.1. Uniform Fields 

Of our 26 GC fields, 5 showed no measurable difference 
in granularity to the simulated random fields (Fig. 1 (left)). 
The stellar distribution of these fields is indistinguishable 
from random in all three wavebands indicating that the 
extincting material has no structure on scales measurable 
at the spatial resolution of our data. The C-C and C-M di- 
agrams of these 5 fields indicate that there is a single value 
of reddening for these stars (consistent with an approxi- 
mate value of Ak = 2.5 for the GC). This single value for 
the reddening agrees with the observed lack of granularity 
across these fields. 

3.2. Granular Fields 

The remaining 21 fields all exhibit some degree of struc- 
ture in their stellar distribution and show a greater range 
of reddening in their stellar population than fields without 
granularity. 

Ten fields show only a small level of granularity, with a 

maximum y value of 0.1 in J, reducing in H and showing 
zero granularity in K^- The C-C and C-M diagrams for 
these ten fields show reddening that is slightly increased 
and more varied than the fields with no granularity (§3.1) 
ranging over ~ 1 magnitude in K^- 

The remaining 11 fields show granularity in all three 
bands. The granularity is highest in J, intermediate in 
H and lowest in Ks\ as expected the longer wavelength 
observations are less affected. Thus the i^g-band images 
will have an apparent measured stellar distribution closer 
to the true underlying distribution. The C-C and C-M 
diagrams for these fields show a wide range of values of 
reddening within individual fields, with the values of red- 
dening 2-3 X that expected for the GC (see Fig. le). 

3.3. Granularity Properties 

The nature of the wavelength dependence of the gran- 
ularity, namely highest y-values in J, intermediate in H 
and low/absent in Ks, indicates it is an effect of intervening 
extincting material. Figl demonstrates this relationship 
between extinction and the granularity in the stellar distri- 
bution. Further, the association of high and varied redden- 
ing in C-M and C-C diagrams with significant measured 
granularity (i.e. high value of y) provides independent ev- 
idence for the relationship between extinction and granu- 
larity. These data indicate that the cause of both observed 
effects is the same; namely, that both arc a result of dust 
and gas structures in the GC. In Fig 3 we plot the range 
of the FWHM of Ak for each field versus their measured 
granularity, these values are given in Table 1 for all fields. 
The extinction values were calculated using the near in- 
frared colour excess (NICE) method of Lada et al. (1994), 
with model giant branch colours supplied by P. Podsiad- 
lowski {priv comm) for comparison. The extinction law of 



Rieke & Lebofsky (1985) was used to calculate extinction 
values from the colour excesses. A detailed spatial analy- 
sis of this relationship, including comparison to reference 
giant branches and analysis of the extinction law towards 
the GC, is presented in the follow up paper (Gosling et 
al.) together with extinction maps of the individual fields. 

The characteristic scale for the size of the extincting re- 
gions is 5" — 15" (Fig. 2), the same scale being measured in 
all three bands in all fields. This angular size corresponds 
to a physical scale of 0.2 pc — 0.6 pc at a GC distance of 
8.5 kpc. The average separation of stars over all fields is 
2.72" in J, 2.29" in H and 1.94" in K^. Fig. 2 shows that 
we can exclude the possibility that the observed granular- 
ity is due to the intrinsic stellar distribution at the 99.7% 
confidence level in all three bands. Note that the y pa- 
rameter in Fig. 2 is truncated below the 3a value of the 
average stellar separation. 

Consideration of the granularity measure in all the fields 
with their positions on the sky showed that the granularity 
does not trace a field's position on the sky with respect to 
Sgr A*; no trend is seen with Galactic I, b, or projected 
radial distance from Sgr A*. In addition, there appears to 
be no relation between the stellar density of a field and its 
position with respect to Sgr A* or its level of granularity. 

We investigated whether there was any clear link be- 
tween the granularity and the extincting material observed 
in surveys at other wavebands. Comparison to ISOGAL 
MIR smvcys (Omont et al 2003; van Loon et al 2003) 
did not show evidence that the granularity traced areas of 
strong MIR emission. A large number of structures ob- 
served in our study seemed to be narrow and elongated in 
appearance, so a comparison of the levels of granularity 
was made with non-thermal filaments and other features 
also observable in the radio (Vallee et al 2004; Yusef-Zadeh 
et al 2004; Mezger et al 1996). No correspondence of ra- 
dio structure with levels of granularity was observed. The 
resolution of published IR extinction studies has been in- 
sufficient to allow comparison to the structures we have 
discovered in our high-resolution (< 0.6") survey. 

4. CONCLUSIONS 

Complex structure in the extincting material towards 
the Galactic Centre has been observed and measured for 
the first time on arcsec scales. The angular scale of the 
structures is 5" — 15" which corresponds to a physical scale 
of 0.2 pc — 0.6 pc at a GC distance of 8.5 kpc. 

Granularity is only apparent in fields in which there is 
high and variable reddening, (derived from the observed 
stellar C-C and C-M diagrams). The granularity is likely 
to be the effect of extincting material obscuring the un- 
derlying stellar distribution. Granularity is higher in J, 
intermediate in H and low/absent in as expected from 
the wavelength dependence of extinction. In 5 of the 26 
fields where a single value of reddening for GC stars can 
account for the observed colours, no granularity is appar- 
ent. 

We remark that the presence of extinction towards the 

GC on far smaller scales than previously observed, means 
that an average extinction correction may in many cases 
not be valid for photometry of individual stars in the nu- 
clear bulge. Schultheis et al (1999) suggested that smaller- 
scale structures in the extinction distribution were respon- 
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Fig. 2. — Statistical measures of the granularity y, in the GC fields. From left to right, these are for J, H, ifg-band observations. The 
rightmost panel is for test fields as described in the text; the darker colours in this panel represent the test fields with higher stellar densities. 
The granularity measured for our GC fields is truncated below S-ct of the average stellar separation. 
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A Ak is the FWHM of the extinction distribution (excluding 
the foreground x population), y is the maximum measured in 
the field. Values of modal Ak are only quoted for fields with 
negligible granularity for which an average Ak is meaningful. 



sible for the observed double-peaks in histograms of stellar 
number versus Ay in the GC. The findings of this Letter 
strongly indicate that extinction may have significantly 
higher values on smaller scales than previously measured 
for the GC. A follow-up paper by Gosling et al., measures, 
maps and analyses the extinction at the GC. 

A very recent paper by Nishiyama et al (2006) has found 
that IR extinction varies from sight-line to sight-line to- 
wards the GC and that the universality of IR extinction 



values are not valid for the GC region. Our paper demon- 
strates that the extinction varies on even finer scales than 
they suggest, in a way which correlates with the granular- 
ity of the stellar distribution in the GC. 

A. J.G. thanks PPARC for a studentship. K.M.B. thanks 
the Royal Society for a University Research Fellowship. 
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under programme ID 071.D-0377(A). 
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Fig. 3. — Relationship between level of granularity (y) in a field 
and the dispersion in the extinction distribution measured in the 
same field. 
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